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ABSTRACT 


This thesis investigates the effects of atmospheric 
turbulence on the spatial properties of the phase of transmitted 
beams. Experimental measurements in the form of a series of 
photographs of interference patterns between two beams were 
made using a modified optical heterodyne system. The two beams 
or the heterodyne system were mixed with a very small, control- 
lable angular difference between tneir directions of avvroacn to 
the detector. The resulting frinse patterns were effected by the 
atmosvhere ses differences between succeeding vhotogranvns snow. 

The first order effect of the atmosphere was found to be a tilt 
relative to the direction of travel. Tilt differences as small as 
5 seconds of arc can be measured using this system with assurance; 
smaller values could be confused with systematic effects. 
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INTRODUCTION 


The advent of the laser has stimuleted interest in conmmuni- 
cation sat ontical freavencies. The potential bandwidtns and ean- 
tenna gains hold great promise, but the degrading effects of the 
‘atmosphere on vropagating optical wavelengsth beams have prevented 
| Immediate widespread usage, Present efforts in the erea of atmos- 
pheric optical communication are mainly directed towards gaining 


an understanding of the atmospheric phenomena and characterizing 





Wthe received fields after propagation through the atmosphere. 





This information ig necessary to build a statistical model of the 
‘channel at optical and near infrared frequencies. <A statistical 
model will allow the use of tools from modern communication theory 
to establish the ultimate limits of optical communication links 
utiliZing atmosvheric chennels. 

Atmosvheric effects on provagated fields at ovtical and near 


infrared frequencies can be broadly put into three catagories: 
lL. Molecular absorption 

Pepa C lem CauLer ne 

3. Turbulence 


| The first two are more or less suggestive end under some condi- 
tions these two effects can be lumped as a path-length dependent 


loss term. This assumption depends strongly on atmospheric condi- 





tions and propagating wavelengths. However, with a clear atmosphere 
‘this model is reasonable at a wavelength of 6328A. This thesis 
deals with the third effect, turbulence. 

Turbulence is the name used to describe the effects of small 
Variations in the local refractive index caused by turbulent mix- 


ing of the atmosphere under wind and temperature gradients. At- 








Oe 
mospheric turbulence results in an inhomogeneous medium, and pro- 
pagation through this medium causes random spatial and temporal 
fluctuations in amplitude and phase of the received field. A 
body of work has been done on the amolitude fluctuations, called 
scintillation, with beams. There has been less empvhasis on the 
phase process due to the difficulty of phase measurements. 

This thesis considers a method of phase measurement using 
an ovtical heterodyne system and evaluation of the experimentél 
measurements. 

Chapter 1 presents a brief survey of the pertinent theory 
and modeling necessery for this thesis. 

Chapter 2 describes the experimental system used to hetero- 
dyne at 6328A and the data analysis system in general. Appendix 
A contains a more detailed explanation of the system, special 
eevices, and electronic circuits. 

“Chapter 3 presents the experimental results. Appendix B con- 
tains more detailed information on the method of data extraction. 


Chapter 4 considers the exverimental results and conclusions. 





CHAPTER 1 


This chavter is concerned wlth some of the theoretical re- 
Sults develoned for beam provagation through the turbulent atmos-~- 
phere with emohasis on those results vertinent to the exverimental 
system constructed, 


1 - 1 Atmosvheric Turbulence 


. 
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The effects of atmospheric turbulence on ovtical wavelengths 
have long been known to astronomers; the twinkling of starlight on 

a windy night 1s a common example. The vossible use of the laser 
for atmospheric communication links has renewed interest in the 

| study of atmospheric turbulence and its effects. 

Turbulence is small scale vatiations of the index of re- 
Mereaction, tyvically on the order of 1 part in 10°, caused by 
temverature discontinuities set uv by turbulent mixing of the at» 

| mosphere. These refractive index changes were first described using 


a Gaussian correlation function which was an analytically convenient 





| method but led to questionable results. A second path taken by 
Tatarski applied the Kolmogorov theory of turbulent flow at very 
large Reynolds numbers to the index be aaeonee This analysis 
leads to a second moment description of the medium, a method that 
is not alone sufficient to describe atmospheric turbulence. Howe 
ever, the fields propagated through the inhomogeneous medium will 
have lognormal statistics; thus a second moment description of the 
medium will be adequate to handle the received field statistics. 
The results of Tatarski's work show that the turbulence is locally 


4 
isotropic and homogeneous and has a structure function of the form 











C (R)p Compete 
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ace ae 
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where 


Pee RS OE a ee I LTE PE ORL AEE. ee Oe 


D(p,H)= (nlz,) - n(x, +e)” 





(2) 


i= En +p | where KR is a vector. <A bar will be 
_ 2 used to denote vectors. 


° 


2 
C08) = constant associated with the magnitude of the 
Diao Lence sab. i 


n(r) = refractive index at r 


Cy = inner scale of turbulences on the order of limmn 
and related to the breakup of_laminer flow. At 
distances greater than fo, n(r) tends to be dis~ 
continuous. 


Lo = outer scale of turbulence; on the order of meters 
and related to the gross layer structure of the 
atmosphere or the distance between obstacles, 


The model, while not as convenient analytically, is senerally 
description 
considered better than the Gaussian correlation function/eand has 


meen the basis for most further work, 





73 =r, 


2 - 2 Propagation through the turbulent atmosvhere 


The available theoretical solutions generally sre for infinite 
plane Wave. sources and assume e "frozen" atmosphere, where the 

ME pulicnce is fixed in space. The frozen-in structure is assumed 

to move across the transmission path due to the comvonent of wind 
normal to the propagation path to account for the temporal behavior. 
The three vasic approaches to the turbulent propegzation problem are. 


using a Rytov aporoximation to the wave equation, solving the wave 











| dw 
a 
equation for the mutual coherence function, and a geometric optics 


3 


fype of solution. The last method, geometrical ovtics, is general- 
ly considered inferior since it does not account for diffraction 
'effects on tong paths with small scale inhomorenetics, 


Meer) Received Field Parameters 


Using the "frozen" atmosphere assumption, the central limit 
maeorem, and one of the solution approaches, it can be shown that 
fa time sample of the received field envelope is a lognormal random 
Variable, and that a time sample of the received field phase is a 
Gaussian random ree blie's There is evidence to expect that the 
Pe on of these statistics is unchanged for beam propagation. 
Another measure of importance is the spatial panesenee of the 
phase of the received field. Tatarski has derived an expression for 
the case of a linearly polarized vlane wave with the atmospheric 
turbulence described as in Eq. 1, and O°) a constant inédevendent 
of a. Then 

e a 
Dp(p) = ack Lp (3) 
where Dol p) is the vhase structure function, k the wave number of 


tne radiation, L the path length, and 
291 p> AL 


ese foe pK AL 
‘Using the expression for Dp(P) in Eq. 3, a coherence diameter for 


Some selected mean-square difference in phase can be calculated. 
hy 


2 


One example quoted in Hoversten from Goldstein is 


_ -2 2 
Beoherent = +*/2 x 10 sale. (4) 
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Fried has used a form of Eq. 3 to show that wave front tilting is 
mace lLarcest contributor to distortion. 

Another parameter is the frequency spectrum of the amplitude 
and phase turbulence processes. Measurements of intensity spec- 
trums have shovm them to be no more than several hundred cycles 

6 
Wide for most ee Measurements for the phase process are 
not as plentiful, but recent measurements have shown that the spec- 


& 
trum is no more than several hundred cycles wide. 


1 = 4 Communication Model 


A suggested model describing the effects of the turbulent 


atmosphere as e@ communication channel,when particle scattering 


mes not important, is a multiplicative channel with additive back- 








q 
Pgeround noise. <A linearly polarized field (chosen for simplicity) 


in the absence of losses, turbulence and background noise would 


be received as 


7 = jwt 
Erect rst) = Real } (F, t)e oo 


When turbulence and background noise are present the received field 


(again in the absence of losses) would be 


jut § (r,t) _—  §@n(F, 6) - 


Beg (tet) = Real } (Ft )e & +: E(rstie e 


(6) 


where 


ob (rst) represents the effect of the turbulent channel 


$(xr,t) = x(t) ig jO( r,t) 
X( 5, mepcesenvs the scintlliation process 


P(r, t) represents the phase distortion process 
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E_(r,t)e 


JO (rs t) 

represents the relevant polarization 
component of the backrround radiations 
& complex Gaussian process with a very 
wide bandwidth on both Hand © | « 


mee= 5 Ootical Heterodyning 


Two waves Simultaneously incident on a surface will cause an 
interference pattern. This can be modeled simply, using plane waves 
of the same frequency and linear polerization. When the two waves 

beat together from slightly different directions, By and Bos 2, 
fringe pattern whose period is related to the difference in approach 
angle is created. If one of the waves has been sent through the 
atmosphere, it will have a phase term proportional to the phase dis- 
tortion process. Thus if complex notation is used, suppressing 
the ce!" time variations for simplicity, the fields can be des~ 


cribed as (see Fig. 1) 
~j8z X(P,t) Jo(r,t) (8) 
e e 


~ JB 2-jBLx 


Ey o(rst) = Ee (9) 
where 
|B — lee: te oc ganee W) 271 —~ -— 
B | = | 8 | = wWre Z Sloe ned Clare 
1 é Be : BY. C r wave bectore’ 


Then on the surface C= eorovmelietie, 1 che total field can be 


(written as 


i 


Feotar (Pst) = Eng (Yet)+ Ey (Fst) 


yd rae) j°(0,t) ~ JB, 
= Eye e “4 Boe (10) 


wt 
—<—, 


where 





Figure il. 
Ooticali Heterodyne 
Geometry 





‘Then the intensity on the surface is 


= Le . > 
I(p,t) = 27 | EtotarZtotar | 


2A( Dig) 2 X(p,t) 


Z 
= 2 13 e “F | 2| +5 2\z| Bol e coal 8a F,t) <a 


(11) 
where is the imvedance of the atmosphere and ua represents the 
'phase angle between the complex amplitudes By and Boe Thus a pore 
tion of the intensity is a term phase modulated by the atmosvheric 
phase distortion. In the experimental section of this thesis a 


System to exploit this effect is described in detail. 





~15~ 
CHAPTER 2 
This chapter describes the system develcved to create a 
Bpeavial enterference pattern between two beams of coherent light, 
one modulated by atmospheric turbulence, 


2- i Systen Layout and Povioment 


The basic arrangement for heterodyning is shown in Fie. 2. 
The source was a Spectra-Physics Model 119 Helium=-Neon gas laser 
operating at a wavelength of 63284. The Model 119 laser has a 
rated continuous output power of O.1 milliwatt although measure- 
ments indicated slightly higher velues during exverimentation. The 
beam from one end wes collimated using a Spectra-Physics Model 331- 


1333 spatial filter-collimator and transmitted over an atmospheric 
| 














path. The down range reflector, shoim in Fig. 3 and Fig. 4, is a 
Sixteen-inch diameter flat mirror mounted in an aluminum collar and 
frame. The collar and frame have precision adjustments througn screw- 
arm systems that allow very accurate pointing of the return bean. 

The beam from tne opposite end of the laser was reflectcd off two 
Fiat first surface mirrors mounted on precision adjustment hoiders 
to reach the combining beam splitter. No lens was used in this non- 
atmospheric path to minimize extraneous lens aberrations. AS this 
bean was close to plane wave conditions by the time it reached the 
combining beam splitter, its curvature was easily matched by the 
telescope output using the telescope eyepiece. 

The receiving telescope is a three-inch reflector with an 


eighteen-inch focal length. The lens from a Praktica FX-3 camera 


(f/2.8, 50 millimeter focal length) was vsed as the eyepiece for 
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Front View of Downrange Reflector and Mount 


Preure 3. 
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several reasons; the 50 millimeter focal length meant the telescope 
output beam size closely matched the local oscillator beam size 

and the focusing ring of the eyepiece allowed adjustments to be 
made with a minimum of mechanical movement. 

The two signals were combined using a fifteen-millimetver 
prismatic beam splitter mounted on a fine adjustment holder. This 
mount was used to manipulate the fringe spacing and orientation 
‘since relatively small changes in angular separation can be made 
most easily at this point. 

All of this equipment, With the excevtion of the large down- 
range reflector, was bolted down on a tyvoeinch thick aluminum op- 
tical bench housed in a penthouse on the roof of Building Beet gl. 'T. 
(See Figs 5 end 6). The downrange reflector was housed ina 
protective box on the other end of the Oot Ofealilding 20, mak- 
ing the atmospheric vath approximately 200 meters round trip. 


fee 2 reattern Detection 


The interference pattern of the two beams was recorded on 
black and white film using a Praktica FX~3 35mm single lens reflex 
camera. The 35mm negatives were developed using Kodak formula D-76 
developer, standard for the Tri-X film used, and were mounted in 
holders to be used as slides for later projection. The negatives 
were used directly for this projection step since the pattern of 
light and dark fringes was the same whether or not the imege inten- 
Sities had been reversed. 

The fringe patterns were examined in detail using a slide 
Beeectox to put the image on a screen. This image was then 


Scanned using a Packard Bell Model 920 television camera. The 
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aes 
Video output of the television camera was displayed on a Tech- 
tronix Model 545 oscilloscope as described in section 2=/!, 


2-3 Exverimental Procedure 


ET TONE IE BS LS 













The alignment betieen the source and receiver was acgompvlished 
by finding the laser (not turned on) image in the receiving tele- 
Scope and down-range mirror system during daylight hours, Thus 
the optical axis between laser, mirror, and telescope wes estab- 
| 
dished. The laser beam could by angular adjustment of the laser 
mount then be put on the mirror at night to have this beam direction 
coincide with the OOCLCal anise 
At this point the beam emerging from the opposite end of the 


Maser was positioned with the small fiat mirrors so that it pointed 


me the center of the beam splitter., All further adjustments were 


made using the beam splitter. The local oscillator signal, making 
&@ right angle change in direction was much more sensitive to the 
movement of the beam svlitter than the beam coming out of the tele- 
scope momie Sstraignt through. Thus the local oscillator beam could 
be moved without disturbing the telescopve output. 

The beams coming out of the beam splitter could be safely ob- 
served by eye after they had passed through a Kodak# 2 Wratten 
filter, Without an eyepiece the beams appeared eas pencil points 
which could be adjusted to coincide. An eyepiece could then be 
used to see the resulting fringe pattern. 

After establishing a fringe pattern, the filter was removed 
and the Praktica FX-3 camera without lens was used to record se- 
quences of three to six consecutive frames of the same fringe pat- 
feetn iec., three to six frames were exposed without altering the 


‘Optical systen. 
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2-4 Data Presentation 


Tne slides used for analysis were scanned using a Packard 
Bell 920 television camera synchronized to an external clock, 
This same clock outout was used to trigger the oscilloscobe, 
The video outout was displayed in two distinctive manners, 


The first method was designed to obtein an overall victure of the 


fringe process across the aperture, the second method to Look at 

Narrow strips of the fringe vorocess across the aperture. 

The first method is shown in block diagram form in Figure 7. 
The horizontal sweep lines, modulated by the alternating dark and 
ight strips of the projected fringe patterns, were displayed one 
on top of the other by setting the time base of the oscilloscope 


| 


to the horizontal sweev time of the television camera. The ver- 









tical deflections of each line represent the tones of grey as 
Scanned. Then a ramp signal was added to the horizontal sweep 
lines to separate them out vertically. This ramn signal was 
synched to the vertical sweep time of the television camera. 

Thus a ouesi-xvelevision picture is made with vertical voltage de- 
Tlections replacing shades of grey. 

The second method used to display the television scan was to 
select out one individual horizontal scan line with the system as 
shown in Figure 8. The "felock delayed by B clock" mode of the 
Techtronix 5/45 oscilloscope was used. The A delayed by B mode is 
@ system where first a ramp Signel set to the B clock time base 
is generated. This ramp signal is compared to a value dialed in 
iy the Peove overator. When the two signals are equal, the A clock 


is armed to take the next trigger it sees and,when triggered, to 
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Figure 7. Full Aperture Display 
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sweep once on its time base. Thus one sweep (A) in a sequence (8B) 
an be selected. By having the B clock synched to the vertical 
myecep Of the television camera and displaying with the A clock 
ynched to the horizontal sweep of the television camera, a select- 


i horizontal sweep of the television camera can be displayed, 





CHAPTER 3 


Chavter 3 describes the experimentally observed results and 
the degscriotions of the phase processes inferred from this data. 


3-1 Beam Descrivtion 


The individual intensities of the local oscillator beam and. 











ithe transmitted beam after propazation through the atmosphere and 





receiver optics are shown in Figures 10 and Il respectively. 

hese photographs were made by simply blenking off the other beam, 
and so revresent typical intensity distributions at the film sur- 
mace for each of mem, The picture of the intensity over the 


‘receiver averture shows the effects of SCiNtitiaciom, 16s, &2 non- 
uniform distribution of intensity with a spatially random pattern. 
The blank spot in the middle and the strip running off to the 

Be ont in Figure 11 are due to the secondary mirror of the Newton- 
‘lan telescove used as the receiver. Although this mirror blanks 
‘off an avpreciavle portion of the vhotogravh, enough of the aver- 
| 

ture is seen to allov analysis. 

When the two beams are present Simultaneously and allowed 


BO interfere on the film surface, the fringe patterns shown in 


Figures 12, 13 and 14 are created. It was found through experi- 
‘mentation that easily discernable changes in tne fringe svacing 


over the averture could not be seen unless the number of fringes 
Was reduced to approximately twenty or less, If the model of Ap- 
pendix A is assumed, then equation (A-12) describes the fringe 


pattern as 
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Figure 10. Local Oscillator Beam Intensity at Detector 
surface 





Figure 11. Transmitted Beam Intensity at Detector Surface 
after traversing Atmospheric Path and Receiving 
Optics. 
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13) = oh Aus” ; thy 2 2X(p) 
of. 


x(p) 
aut| Ua) | Ue col Be -p()-a3} (Aw?) 


U. revresents the amplitude of the envelope of the field 
sent through the atmosphere 


Us represents the amplitude of ‘tthe envelope of the local 


oscillator field where a is the phase angle between a 
5 and. Bo 





N represents the average collected background rediation 
Intensity 


i represents the losses in sending a beam through the 
atmospheric channel and in collection 


B. represents the direction cosine between the two waves 
at the detector surface (see Fig. 1) 


y represents the impedence of the aetmosnhere 


UW represents the magnification of the telescope and eye-~ 
picce 
axle) 
represents the scintillation process due to turbdulence 








(5) represents the phase distortion process due to turbulence 
This intensity expression reflects the fact that this pattern is 

assuned to be a time sample of the fields. It contains only spae- 
tially varying terms and implicitly assumes that this pattern is 

Maken for a particular time t. At some later time a photograph 

would have another time evoltion of the sample functions X(p) 


nd D(p). 


As shown in Apvendix B, the geometry can be selected so that 


the fringes occur along one of the axes in the p plane. With this 


meometry the direction cosine, By? can be related to the difference 


in approach angeles, f » of the two beams so that the intensity be- 
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2e sequence A Fringe Patterns. These consecu 

_—* Bere ec abhs will be analyzed in the followin 
figures: Figs. 


22 
25 





= Or 





Figure 13. Sequence B Fringe 
Patterns. These three con- 
secutive photographs will be 
analyzed in the following 
figures; 

Figs. 16 
a? 


23 
25 
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Figure 14. Sequence C Fringe Patterns. These consecutive 
photographs will be analyzed in the following 
figures: Figs. 17 
20 


24 
er 








comes 


where 
| 


G@ = the phase angle between FE, and E 
) sin? = 18 ,| 718 | ~ Q 

_ ant 

It can be seen that the fringes are caused by the phase modulated 
Binvsoid riding on the intensities of the individual sisnals. 


murcvher, it can be seen that if the period of the sinusoid, #, 


is small compared to the distence over which significant changes 






occur in the values of Dl), the ohase distortion, individual 


Tringe lines will not experience much perturbation. Since a 









maximum of twenty lines distributed over the three inches of 
perture show avpreciable perturbations of the fringe lines, a 
minimum significant distance in D(5) can be estimated at (3/20) 


maches = 4 millimeters. This rough estimate of significant di-~ 


general experience from longer vaths. 

, The amplitude scintillation in the three sets of fringe pat- 
erns does not appear to be severe, but this is due to a selection 
meocess for just this quality. 

Another notable feature of these slides is the spatially vary- 
Ing fringe pattern across the Sener The fringe spacing in 


Figs, 12 and 13 are narrower on the left side of the aperture. 


Especially in the photogranvhs of Fig. 13, one can see that the 
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fringe spacing not only varies spatially, but also chances radi- 
cally from one photogravh to the next. The photogravhs in Fig. 14 
show a rotation of Lines between photosravhs. The photographs of 
Bie. 12 were taken just prior to those in Fig. 13, while those in 


Fig. 14 were taken on another night without any appreciable modi- 
















fication of the system having taken vlace, 


1 


A photograph of the system through a non-atmospvheric path 





could not be obtained unless the system as installed was changed 
significantly. The divergence of the wave emerging from the laser 
or laser-collimator when beamed into the receiver telescope was 
such that its focusing point was moved too far out to use the ca-~ 
mera eyepiece as mounted. This camera eyepiece was needed to match 
curvatures with the LOwave. Other eyepieces required optics in 
theibO path. When the divergence was reduced so that the system 
could beevuced “as is", the aperture was not fully illuminated; and 
very little of the averture could be seen around the secondary mir- 
mor. Thus system effects can not be directly taken out. However, 
the availability of consecutive photographs allows analysis of the 
changes in successive photographs. The ovtical system is the same 


between vhotogravhs exceot for advancing the film in the camera, 


= + a Pk, IO ee! RE Co 


In an effort to obtain a feel for the vrocesses that occur 


Over the entire averture, the slides were projected on a screen 
























a Bl pn 
and scanned with a television camera. The output of the television 
Camere was displeayed on an Oscilloscove using the first of the 


methods described in section 2-4, This method is essentially 


Se ac 
the construction of a television picture but with the shades of 
grey reolaced by vertical deflections on the individual scan Lines. 
This method gives the appearance of a three dimensional mav of the 
c 
intensity over the aperture and the results are shown in Figs, a 
16 a 17. 


The main interest in disvlaying in this manner ig to ret 


mome feel for 


(1) Phase changes along the fringes 
(2) Phase changes across the fringes 
(3) Seintillation 

(4) Systematic effects 


The deev troughs on the edges of the photogravhs in Figs. 


= 16 and 17 are the locations of bleck non-reflecting tape put 


{| 
on the vrojection screen to sid in centering the picture. In e- 


very photograph, the fringe svacing increases towards the right. 
mois is considered a system constant and as such will be dis-~ 
counted in further analysis. 
| The fringe lines tend to curl towards the lower right-hand 
corner. This curl is stronger in the lower half of the aperture 
mad can be seen particularly well in photograph of Fig. 17. This 
is an example of the phase distortion along the fUringe lines. A 
Islightly different effect is the amplitude modulation along the 
fringe lines caused by phase distortion. 

The scintillation can be seen best on the peak lines as on 
amplitude modulation. Where the modulation produced by the fringes 


is shallow, generally on the right hand sides of the pnotogranohs, 
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Figure 15. Sequence A Full Aperture Display of the 
Fringe Pattern of Photograph#l anda2 
by Method of Section 3-2. 
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Figure 16. Sequence B Full Averture Display of the 
Fringe Pattern of Photographs 1 by Method 
of Section 3-2. 





Figure 17. Sequence C Full Aperture Display of the 
Fringe Pattern of Photograph #1 by Method 


of Section 3-2 
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scintillation tends to swamp out the detail of the Mattern. 






The analysis of the phase vrocess across the fringe grain 
mill be done in detail in the next two sections. It is worth- 
Nile noting at this one that the frinse pattern is generally 
Meeto.e and enough definition exists to justify a Look I Cevai 


5. 53 Phase Analysis 


Using equation (B-2) as the model for this fringe line sys- 


prem, the intensity is 


SS 1 c 2 2 gis. 
I(p) 344 up| + N+ ea 
(p) 
2TtyYxX 
+ alu | Poot -fergrath 8 
if a narrow strip in the y dimension is examined, narrow enough 
SO that the processes do not change with y, the intensity ex- 


pression becomes 


| I(x,y,) = Atel + | + <e in tule 


P out, 2, on ¥0) poor Bul ~ bix.yg)-<f (13) 


Assuming that the spatial sinusoid is visible despite the random 
amplitude modulation from scintillation, the distance between suc- 
cessive peaks or valleys can be measured. This indicates the sin- 
usoid argument going through 2m radians or one veriod. Thus, 
choosing the peak lines as the point to measure period length on, 
we know over one period 

see -P(%s¥Q)) - (2m | boty) = 2m (14) 
where | 

x and x' are the points where two adjacent maximum valves 

of the sinusoid are Located. In the remainder 


of this thesis x and x' will indicate adjacent 
maxima. 





<BR 


blince x and x" can be measured, we get 


) pe 
Ply) -Uxtyyg) = enlE(x-x") = 1] (15) 
Thus over the distance (x-x') we get a measure of how the phase 


Jistortion process acted. If we continue this analytic procedure 






across the fringe pattern, we have sone notion of how (x,y) 
mianses with x, 

lt has been shown by Fried that geometrical Boe oe 
xan be placed on the form of (x, vie ). If a Taylor type expan- 
ion is made of D(xsy0), we get 


D(x, 9) = A+ Bx + Ox% + caer (16) 


rutting this expression for D(x4¥6) into Eq. 13 gives 


1 Z a Zee 2 2X( XY.) 
I(x,¥,) a aA ieee \U, e 


KGa ) 
° cos[ STAs - C - A = BX - Cx”, os \ 





+ euL| | Z| e 


Oh 





jere it can be Seen readily that 


A represents a phase origin relative to the position 
coordinates x 


the direction of travel 


2 : P 
Cx” . . e represents a spherical phase change, representing 


@ svherical wavefront on the received wave 


| Bx represents a linear phase change, a tilt relative to 
r the coefficients give a good measure of the importance of each 


erm; the number of terms needed to adequately represent the phase 





mrocess indicate in some Sense how simple the vhase chenges are. 
There are two problems involved with applying this analysis 
lirectly to the fringe patterns made for this thesis. The first 


| 


is that no non-atmosoheric vattern exists to take out the purely 
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systematic effects. The only method of removing these systematic 


effects is to examine consecutive patterns and remove the constant 
Bi fects: this leads to the second problem. The only distinctive 
feature that can be used for locating the same narrow strip in 
Bonsecutive patterns is the secondary mirror. Thus without taking 
a Por several entire slides, which would take a prohibitive 

| 


Mounmeeor time, only limited strips can be handled in each pattern 





















Or comparison with the next. Fortunately, the regions just above 
and below the secondary mirror are not unique compared to the 
remainder of the. fringe patterns in the rest of the aperture, as 
san be seen by Bee he Figs. 12, 13 and 14. Thus these regions 


kre & reasonable choice for analysis. 


} - 4 Narrow Strip Anelysis 


To examine how the fringe vatterns made for this thesis 

Be nzed in phase, narrow strips above and below the secondary mire 
‘or were measured in detail. The slides containing the fringe pate- 
lerns were projected on a screen, then scanned using the televi- 
Boy camera, The fringes were aligned so that the horizontal 

kcan of the television camera was across the grain of the fringe 
attern. Selected horizontal lines from the television camera 

ere displayed and the locations of the peaks of the fringe lines 
= measured for five to fifteen successive horizontal lines. 

‘his corresponds to a striv one to three millimeters wide over 

Ihe objective aperture. Typical displays of these horizontal lines 
xe shown in Figs. 18, 19 and 20, where the top line corresponcs 


 @ scan line in the top half of the fringe pattern; the mirc le 


4ne to a scan line in the middle, and the bottom line to a scan 











Figure 18. Sequence A 
Typical Horizontal Lines 
from Television Camera 
for Photograph ql 
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Figure 19. Sequence B 
Typical Horizontal Lines 
from Television Camera 
for Photograph#l 


Figure 20. Sequence C 
Typical Horizontal Lines 
from Television Camera 
for Photogranh #1 








alae 


mine in the bottom half of the fringe pattern. The relatively 



















mati modulation in the middle of the middle lines is due to the 
secondary mirror. The relatively deep excursion of the scan lines 
iv the edges is caused by black non-reflecting tape vlaced on the 
reen to provide a means of centering the fringe pattern on the 
elevision monitor. The spatial scale of the oscilloscope display 
an be related directly to the spatial scale of the averture. The 
onversion is 0.40 inches of averture is equal to 1 centimeter of 
} oscilloscope Gisplay. This conversion will be applied to all fu- 
mre figures that deal With the averture. 

Tyvical unavereged strips are shown in Fig. 21, where the 
leak lines are draim in. The locations of these lines are simvly 
weraged slong the fringes. The distance between the averaged lo- 
ations of the peak lines across the aperture is plotted in Figs. 
2, 23 and 24, From Eq. 15 this distance between peak see can be 


felated to the ohase distortion and fringe period as 


(x-x') = ae [O(x,y,) - Plxtyy,) + 2m] (18 ) 
fhus measuring (x-x‘') gives a number proportional to p over 

he distance (x-x'). Further geometric interpretations about the 
orm of c can be made as successive (x-x') values are measured and 
ted across the aperture. <A direct measure of the spatial 
eequency Eig not available, however this term can be assumed 
onstant across the aperture. Then any changes in the (x-x') 
cross the aperture can be attributed to changes in P (x59) and 


| 


ny changes in (x=-x') between photographs can also be attributed 


© the change in D (x49) with time. 


This method of vlotting the average distance between peaks 
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aA 

Rives a visual disvlay of the general type curvature of tke phase 
process across the aperture and seems to show @ general quadratic 
Torm indicative of spherical distortion of the incominse wave, 

Another method of plotting is to use the averazed location of 
pne peak lines on the abcissa and to plot the 2 radian change of 
| 
the Sinusoid argument between successive peaks on the ordinate, 
‘This is shorm in Figs. 25, 26 and 27. These grephs have a linear 
mean square estimate line loereamter wererence. The first aspect 
to note is the dominance of the linear term in every piot. MThe 
a have some higher order deformations, but the linear term is 
‘obviously most important. The second aspect to note is the chenge 
in slope between photographs in Fig. 26. This means 2 significant 
change in tilt of the incoming wave has taken place. We can 


equate the argument of the sinusoid in Eq. 13 to the estimete lines 





Min Figs. 25, 26 and 27 as 





amp ox 
[7 = Dlx,¥,)-0] = yx + (19) 
where 
= atl the linear estimate line slope measured for 
a ” 2 each photograph 
oe 


May the covariance between the averaged location of the 
peak lines, x, and the total phase change, y, at 
moat x 


O the variance of the averaged location of the peak lines 
m the mean of the averaged locations of the peak lines 
mM the mean of the total vhase change across the aperture 
ae 


Ko = My -——5 mi the linear estimate intercept for a 
| os photograph 


‘The change of this argument between successive vhotogravhs is 


P(x,¥.) -- Pixiy.) = (k,-k", )x - (ko-k',) (20) 
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Aperture Sequence A 
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Aperture Sequence B 










NARROW STR? ABove 
27 ECONDARY MiRoel 





- Y 
= 


tS 
= 3 





TOTAL PHASE CHANGE ACROSS 
AFERTORE CRADIANS) 
3 


Mi eee ee eee 
Oo of O®8 12 £6 #420 2s 2.8 
DISTANCE ACReSS APERTUCE Cincres) 


© DATA POINTS FRoM Fiesr RismaArn 
QO DAG Posts F2oM Ssconp Ps DSCAPA 
A, DATA POINTS F204 Tien Prsracart 









NAZRow STRO Beran 
SECONDARY MiReeR 


ga 8 


ASE CHEAIGE ACKOSS 
CIF DUIS) 
S 


APERT ORE 
aes 


HH 


TOTAL PF 
at 


I ac A 


fae ee eee 
6.0 04. 0% ae i, 20 2.4 2.8 
DISTARCS Ackoss APERTVEE CINcHES) 


Figure 26. Total Phase Change Across the Receiving 
Aperture Sequence 8 
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S he _ 
where the primes refer to the succecding photograph. The (k,-k° 
LS not meanineful since many steps in the analysis introduce arbi- 
brary references in phase. However, the change in slope is easily 
measured. These changes in Ky ace shown in Table 1 and the eauiva- 


ment change in tilt of the incoming wave is also listed, 





| Ax(radians/inch) Tilt chanse (in 
seconde O2 ate) 
Above secondary mirror Nee See 
igure 25 
Below secondary mirror Za ON as 
Above secondary mirror 4.251/.75n 10,907 lees 
igure 26 | 
Below secondary mirror 8.751 Cia 
| Above secondary mirror (Sn » 386 
Bigure 27 
» Below secondary mirror ee 13925 
Table 1 Measured Changes in Tilt in Three Sets of Consecutive 


Pieesographs (from Figs. 12, 23 ana 14) 


As is discussed in Appendix A, tilt changes less than 15 










microradians (3.2 seconds of arc) will be indistinguishable from 
possible mirror tilt caused by vibrations. This limit is exceeded 
Only in Fig. 26. If the photogranhs used to generate Fis. 26 are 
Peexamined (see Fig. 13), it can be clearly seen that the fringe 
macing changes dramatically from the first to second photograph, 
Brfeot Mec seen in Figs. 12 and 14. 

fn an effort to determine the higher order deformations, a 
fOs—-log plot was constructed using the average location of peak 
mines as the abscissa and the vohase changes,with the linear esti- 


Beve subtracted out, as the ordinate. However, oniy systematic 


meeccts are visible, 
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This chapter summarizes the results presented in Chaver 3, 
and draws some conclusions about the general procedure. Some suge 
fg estions for imorovement of this system are also made, 


| 


Be ~- i Conclusions 


Using the methods outlined in Chapter 3, it is clearly pos- 
Bite to measure phase in a narrow striv across the objective 
Bacerture, Using the resulting phase plots, it 1s possible to com- 
rte a change in tilt between successive photozrapns. in this ex- 
periment, the linear terms and systematic effects masked possible 
nigher order deformations caused by atmosvheric turbulence, but 
With more care on a longer vath these effects shovla be measurable. 

The details of the analysis show that the minimum tilt changes, 
sevsrate from mirror induced tilt, measurable by the present sys- 
tem was avovroximately 5 seconds of arc. It is interesting to note 
that the measured tilt ranged as high as 22.5 seconds of arc, a 
large value for a three inch averture and a short path. 

When the fringe pattern produced by the two beams is looked 
et through a Wratten filter, the system dynamics can be seen graph- 
Bealty. The fringe lines waver along their length and, a vs their 

Spacing os wider, a bowing travelling across the pattern is of- 

men seen. A third effect is a shimmering of the entire pattern, 
Possibly indicating a time change of phaseorigin suggested in 
eo" 3-3. 
. The system of oscilloscope display and visual reading is 


Srearly a tedious method of data taking. The amount of data 


Chat can be taken in this way is small, and is not compatible 








ee a 


ats 


with machine computation. Judgement can be exercised to pick 
Out the peak lines from scintillation when this display method is 


used, but the peaks and valleys of the fringe vattern are the only 























Darts that can be used for visual data taking. 

The lack of a system-only series of photogravhs was a hin- 
drance to the enalysis.— Some further feel for the motions of the 
down-range reflector would have been invaluable in interpvretating 
the results. 


4 vugsestions for Improvement 


The first suggestion would be to set up a path in the la- 
boratory long enovgh so that the divergence of the wave entering 
Dhe telescope is not too fast. At the same time, the beam should 
ijjluminate a large amount of the averture. Thus, using the same 
System as one would use when taking atmospheric data, a recording 
bf the system alone could be made. This would be most helpful in 
analyZing single photogravhs, and sequences would not be necessary. 
A second suggestion is to convert the video output of the 
Welevision camera to digital data for later machine analysis. The 
Dasic equipment is available, and only a gating circult and san- 
ling ae added, The master clock that generates 

me pulse train used to synchronize the television camera runs at 
mO,250 Hz;. the first frequency used in the television synchroni- 
Wation is at 15750 Hz. The sampling pulses could be selected to 
mori this master clock with some flexibility as to the desired 


meequency, and synchronization with the television camera could 













tailed analysis over the entire averture. 

As longer paths are used, the problem of reflector induced 
tilt would be lessened. The increase in beam diameter would not 
be too large as path length is extended, but beam motion due to 
reflector tilts would be more apparent. The motion of the beam 
could be observed and data could be taken only when the beam is 


Steady on the target. Thus mirror induced tilt would be minor. 
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APPENDIX A 


This appendix contains an analysis of the heterodyning sys- 



















bemein more detail than in Chapter 2, a fuller exvlanation of the 
data analysis system, and electronic circuit diagrams. 


g « s 
= | Heterodyne System Analysis 


The scheme outlined in Figure 2 for creating a spatial in- 
werference pattern proportional to the vnase distortion process 
caused by atmospheric turbulence can be analyzed using simple 
models. 

The front end of the laser (typically the output of the 
Spectra-Physics Model 119 laser would be 0.1 milliwatt at the 
front end and 10 microwatts at the back end) is collimated and 
transmitted through the turbulent atmosphere and then reflected 
back over the same path by the down range reflector. After cole 
lection by the receiving telescope tnis signal is mixed with a 
ocel oscillator signal taken from the back end of the laser. 
These two waves can be avproximately described before transmission 


as plane waves of finite spatial extent. As the finite spatial 


experiment, it is not explicitly included. Thus in complex nota- 
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The atmospheric effects on the propagating beam include dif- 














fraction, scattering, and absorption losses; atmospheric turbulence 
modulation; beam spreading due to imverfect reflector-mirror flat- 
ness; phase modulations due to mirror motions caused by vibra~ 
tion, Collection of the returned beam will introduce further 
sosses due to the finite collecting averture which is smaller than 
the beam size. Background radiation will also be collected in the 


receiving telescove, Tne various effects will be described as fole 
e represents the turbulent channel 


y(r,t) = Beret) i jQr,t) 
X(r,t) represents the amplitude scintillation process 
Or, +t) represents the vhase distortion process 
L(r,t) revresents logses due to 
ieee err raction 
2. Particle scatterings 
3. Molecular absorption 
4, Mirror imperfections 
5. Collection losses 
) 


HO (rete - represents the wide bandwidth Gaussian noise 
collected by the receiver 


jQ FROEGsents chase modulations due to mirror movements 


Ising the channel model of Chapter 1, the transmitted field will be 


received as ite an — 
yr jut Yio, oy ieee) 
e. e€ e 
Q(x Yet) jwt 
ey (r,t) e e (A-3) 
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ine collecting telescope will minify the portion of this re=- 





| we BL baw 
celved beam within its aperture. As long as one looks "close 
enough" to the eyepiece, the minified image is a clean represen~ 


tation of illumination in the object space. (Going further away 













gets into the Fresnel Zone of the secondary mixror and aperture 
edge effects.) The collected field is mixed with eae local oscil-~ 
lator field by a bean splitter snd projected on a plane screen for 
detection. The geometry of this plane screen (Fig. 1) is selected 
so that Zz=0 on the surfaces the transmitted wave approaches the 
surface normally and the local oscillator signal approaches the 
surface obliquely, and thus has a direction cosine on the surface. 


The fields using the reometry of Figure 1 are 


~§B ep jut 


Erol Pst) =e Use € (A-4) 
_ " y¥( pst) jut 
Errang(?s ©) = L(p,t)uu,¢ € 
_ j€, (pt) jsut 

+ HE (p,t)e e (A-5) 

Esotal = ErolPst) + Enrang! Pst) (A~6) 
where 

— = the two-dimensional vector in the plane of the 


detector surface 
the collecting telescope-eyepiece magnification 


E 
i 


mhe detector, in this case film, will measure the intensity which 
(E , 

ane (Et otal? total 

( 0 5 1 ) — 9 > > / ae 


time averaging over the period of exvosure as shown in Fig A#1. 


is the veriod of the shutter ovening and is on the order of 2-10 
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4seconds for che 35mm camera. The inpul sntensity in Vee 
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FIGURE At 
. BLACK BOX MODEL OF DETE CTOR 


ter period compared vo ae 


Due to the shortness OL mule shut 


atmospheric process® atl he terms except chose containing noise 
factors are constant over the sanple oer os. i Oy the noise terms 
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will vary; thus the antegratving 
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The integral of the first term, the square of the envelope 














of a comvlex Gaussian vrocess, is a relatively long-term integra-~ 
tion of a very wide bandwidth vrocesss; thus the time integration 
can be replaced with the ensenble average 
mn 2 a |S Pee 
= x Jence,e) dt E| 2 lea | = NN 
where nc is provortional to the average power in the relevant 
background rediation. 

Similarly, the time integrals of the terms containing both 
E(p,t) and © (65%) may be revlaccd with ensemble averages since 
both EL (Pst) and O(p,t) are very wide bandwidth vrocesses com- 


Bared to the integration time T. Thus, for example 
{ T ] 
=a \ ape UE, Ft) Cosl Be F- en(fO]olt = F | 4 Ue ee cos | BoP - onthe] | 


)Further, the two processes HP st) and © (p,t) are independent, 


SO that 


E[ 2). Extocostfe?- eueal]= EL 2*¥- Eattd | EL esl Brf- Ongsl 


Also,®_(p,t) 46 reasonably modeled as approximately uniformly 
distributed between O and 27 at any time during the veriod O to T. 
Thus E | cosl8- 5 - Q.(o,)]| tends towards zero and the time inte-~ 


Beal is aporoximately zero over the entire surface, Similarly, the 


P ee XCe, - 
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term tends towards zero. 

This integration model of the detector imolies that the de- 
tector records energy. For this model, since the vrocesses are 


Slow compared to the shutter time, a typical intensity, I(9), has 



















pbeen used with the understanding that multiplication by T, the 
Integration time constant, will yield the prover energy terms. 

ii@e resulting output from this integration model is thus 

. 1 2 DD 22) 
i)! = a7 ; luo] + ub (6) |v, | e + N 
a: x(n 2 
+ 2uL(pe) | Us|] Bale cos Bp - ){b)= Ye) if (A-11) 

The loss term L(p) is reasonably modeled as a uniform quan- 
Tity over the strong portion of the transmitted beam since scat- 
tering and molecular absorvtion are essentially non-preferential 
and are less aovarent in the strong portion of the beam. As it 
is this portion of the beam that will be collected, the term L(p) 
will be replaced with L. 
The term w(p) represents the phase modulation due to mirror 
vibrations further acted on by the turbulent channel. Due to 
Mirror mount massiveness, the value of w(p,t) is constant with 
Mtime over the shuttering period. The svatial variation of Wp) 
Mmeerucacal, however, Since the mirror movements that can occur 
muse a tilt on the reflected wave. This tilt is of the same form 
miat one expects as the first order effect of atmosphere dis- 
tortion. 
There are arguments that seem to indicate that while this 


milt caused by the mirror is a limitation on the system, the mea- 
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surements of the atmospheric effects can still be made. The first 
is that tyvical tilt changes betvieen photozravhs are on the order 
of 15 microradians. If this tilt change were due only to mirror 
movements, this would mean a beam movement of 0.05 inch at the 
Ptelescove., The overlap of the returned beam over the receiving 


_ aoproximately 0.5 inches, meaning that this movement could 











mnot be seen. However, more Significant tilt changes would be 100 
mpius microradians, which translates to a 0.5 inch travel at the 


jtelescope. Thus this movement would be seen and neglected. Thus 


bath do not swamp out larger values of atmospheric tilt. There- 


more, v(o) will be lumped in Op), and. the intensity can be written 


GG) = zy 0, oe ul 2 axe) 


X(p) 


+ zu u,| | Ure cos[ 8). Op = b()-«)} (A-12) 
where the intensity contains a term phase modulated by the atmos- 
Dheric phase distortion. 


~_2 Electronics 


The data analysis system was built around the Packard Bell 
920 television camera and is shown in block diagram form in Fig. 
m-2. <As produced, the camera has a 2:1 random interlace which 


Mas altered to @ single vertical scan to allow accurate reading 


. 

L. the fringe pattern. The sweep of the television camera was 
lade into 256 horizontal lines each vertical sweep, with synchron- 
iZation between the horizontal and vertical rasters. The circuits 

used for this system were built around RTL digital integrated= 


lo 
mircuits. 
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Figure A-2 Data Analysis System Block Diagram 
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4 master clock shown in Fig. A-3 was used to generate the 
pulse train for the system. 
A series of counters and pulse shavers were used to generate 
synchronization pulse trains from the master clock signal. In 


each case a synchronous counter was used to minimize countdown 


















‘delays innerent in sequential counters. Figs. A-4, A=-5 and A-6 
fare circuit diagrams for these counters and pulse shapers. In 
eacn of the circuits the counting is done by the JK flip-flops 
gated on and off by the logic. The shapers were pulse generators 
With a variable resistor used to set the output pulse width. 

The synchronization pulses were put into Schmi¢t trigger 
circuits to provide some vower gain and also to buffer the count- 
Erecsrcuits, 

Using the system outlined above, the television camera scan 
was made to cover the 256 horizontal lines in each vertical sweend 
with no interlace. The starting point for each vertical sweenp 
was the same so that the video output could be disvlayed on an 
Becamroscope, using the synchronization pulses as triggers, with 
Hprecision and stability. 

Before disvlay, the video output of the television camera 
e's put into a low vass filter with a 3 db point at 4 MHz. The 
mlow pass filtering was done to remove a large amount of high 


frequency hesh. Fig. A-7 shows the improvement of the display pro- 


Mduced by this filtering. 
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Figure A-3 Master Clock and Pulse Generator for 
Synchronization System 
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Figure A- Seven Counter and Pulse Generator 
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Figure A-G Four Counter and Pulse Generator 
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Figure A-7. Improvement in Picture Quality from Insertion 
of Low Pass Filter at the Television Video 
Output. Uppver Picture is without filtering; 
Lower Picture is with filtering. 
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APPENDIX B 


nape @ 8805 A Om ee 


This appendix contains a discussion of the data analysis sys- 
tem limitations, 


B= i Fringe Pattern Model 


As shown in Appendix A, the intensity for vilane wave models 
fis = 

. rn Z 2 2a Ze 2x®) 
ri)» af wal 2 8 9 wat 


+ 2uL| Uy | UI e cos[B y 0 ~- p(F)-c)f (A~12) 













The fringes result from the spatial sinusoid, which is proportional 
to the product of the two signals, The fringes ride a bias which 
is the sum of the intensities of the two signals. The terms re- 
presenting atmosoheric turbulence effects can be thought of as 
minified revresentations of the processes at the objective of the 
receiving telescove. Thus 

igen) = X(p'/u) 


Op) =O(o'/u) 


Pp is the two dimensional vector at the detector surface 
p! is the two dimensional vector at the objective plane 
oriented so that the geometry, when provagated through 
the receiver optics, matches that of 6. 


u reoresents the magnification of the telescove~eyepiece 
combination 


Further, aligning the geometry of this model so that the fringe 
lines coincide with one of the axes in the © viane Simplifies the 


Motation. With this assumption, the intensity is 
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Peep ) 


where 
a [6] (6, Ale) = |8|sin} = sinks 2uL/d 


q represents the angle between the two interfering waves 
and is set by the exverimenter so that «1 


By booking at a very narrow strip in the y dimension so that vari- 


-ations in y do not have a chance to occur, the intensity model can 
be written as 

2 2) 2 PX%S 
I(p) = s+ \y 1 +N +uL |Up| e 
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Xl ey ) 
' f> O a ankx rT «xs «o 
-- 2uL|U,} |U,lc cos[ x ~ P(%4¥0) aif (B-1) 
The problem of seeing the fringe vattern is related to the 





depth of modulation, which can be shown in this model to depend on 
the visibility ratio V, defined as 


vy = amvlitude of sijusoid 
amplitude of bias 


X(p) 


a uh} |22] » 
2 - yal gery) (Be=2) 
i | Y ee lu,| e 
hor this 6briment, i ae henligible compared to the other terms 


in the denominator end will be neglected. The other terms in the 
denominator represent the local oscillator and received beam in- 
tensities. The shedow pattern of these intensities can be seen 
during experimentation. From observation, the shadow pattern of 
the received field is brighter then the local oscillator by an 


estimated factor of three to four. On @ longer path, losses will 


surely be more severe and this ratio would be reversed, If Eq. Be-2 
| 
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is manivulated slightly, and the information on the ratio of 


the two bean intensities incorporated, the visibility ratio be- 


comes 
V = 2 (us [Ua] / Bol Ve 
2 
ieee L “\u,| 2xX(p) 
sh 
yy | 
= 2 KeXlP) ye Xt P) 

Tn eee) | ee (B-3) 


It can be seen that the scintillation process is the fundamental 
‘limitation on modulation devth for this path. This is plotted in 
rig. Bel e 
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Figure Bel Plot of Visibility Ratio with changes in 
the value of scintillation 


merurcther limitation peculiar to this particular data 
analysis system was the resolution available using an oscillo- 


scope display for detailed analysis. The scope face wes ten 





centimeters across and distinctive vortions of the waveforn 


could be located accurately to one-half millimeter as long as 
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ithe surrounding clutter was not too bad. However, at least a 


[three to five millimeter separation was required before 
Weecible to reliably distinguish peaks from clutter for 
cal depths of modulation produced py the scanned fringe 
‘This means that 20 to 33 lines maximum could be handled 


scope. It wes found that (see section 3-1) the effects 


Lewes 


the typi- 


On tie 


atmospheric turbulence on the fringe perturbations did not become 


easily visible until the number of fringes across the aperture 


became twenty or less, Thus the use of the scope was not a 


Serious limitation for this experiment. 
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